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Abstract. We discuss the structure of decoherence-free subsystems for a bosonic channel affected 
by collective depolarization. A single use of the channel is defined as a transmission of a pair of 
bosonic modes. Collective depolarization consists in a random linear U(2) transformation of the 
respective mode operators, which is assumed to be identical for N consecutive uses of the channel. 
We derive a recursion formula that characterizes the dimensionality of available decoherence-free 
subsystems in such a setting. 



1. Introduction 

Decoherence induced by uncontrolled interactions with the environment is a 
major obstacle when implementing protocols for quantum information processing 
in real physical systems. The effects of such interactions can however be reduced 
by using particular quantum states that are robust against specific decoherence 
mechanisms. An interesting and physically relevant scenario of decoherence is 
when the physical system exhibits certain symmetries in interactions with the en- 
vironment. Such symmetries imply the existence of whole subspaces that remain 
completely unaffected by decoherence and can therefore be used for faultless quan- 
tum information processing. This observation, which can be made for a number of 
interaction schemes from various perspectives ^ El [3] , has led to the develop- 
ment of a general theory of decoherence-free subspaces and subsystems, reviewed 
recently in [5]. 

The theory of decoherence-free subspaces and subsystems applies to a vari- 
ety of physical systems. One such system is light travelling through an optical 
fiber, which usually exhibits random birefringence caused by fluctuating environ- 
mental conditions, such as temperature and mechanical strain. Mathematically, 
this system can be modelled as a bosonic communication channel. A single use 
of this channel is defined as a transmission of two bosonic modes corresponding 
to orthogonal polarizations. Birefringence is described as a random U(2) trans- 
formation between the operators of these two orthogonally polarized modes. In a 
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realistic situation this transformation remains practically constant over many con- 
secutive uses of the channel, for time intervals comparable at least with the round 
trip time in a fiber [Jj. The presence of such a symmetry leads to the existence of 
non-trivial decoherence-free subsystems. Their structure has been discussed in our 
recent paper [H]. In this contribution, we will concentrate on the mathematical de- 
tails of the derivation briefly sketched in the previous paper, and will illustrate the 
presentation with diagrams that provide additional insights into the calculations. 

The decoherence-free subsystems for a collectively depolarizing bosonic channel 
have recently attracted considerable attention in the context of quantum reference 
frames and various scenarios of quantum and classical communication over depo- 
larizing channels 0. This has produced theoretical proposals for robust quantum 
key distribution ^Qj ^2 El , alignment-free tests of Bell's inequalities ^H] , and 
entanglement-enhanced classical communication ^1] , with some of these ideas im- 
plemented in proof-of-principle quantum optical experiments |15| I16| I17j . From 
the mathematical point of view, an interesting aspect of the collectively depolar- 
izing bosonic channel is that in contrast to particle-based scenarios in which every 
elementary system (i.e. a particle) has a finite number of states, the elementary 
system in our case (i.e. a pair of bosonic modes) is infinitely dimensional and 
its transformation under the interaction with the environment is described by a 
reducible representation of the group U(2). 

This paper is organized as follows. We start with a discussion of the polar- 
ization transformation for a single use of the bosonic channel in Sec. 2. We then 
analyze in Sec. 3 the emergence of decoherence-free subsystems for multiple uses of 
the channel. In Sec. 4 we derive a recursion formula for dimensions of decoherence- 
free subsystems. Finally, Sec. 5 concludes the paper. 



2. Polarization transformation 

We assume that the quantum system transmitted in a single use of the channel 
is composed of a pair of bosonic modes corresponding to orthogonal polarizations, 
with the respective creation operators labelled by a} H and Oy. Random depolar- 
ization consists in a linear unitary transformation of the modes given by: 

(^)^n(^) (i) 



where the matrix 



»h £ "C i 



is an arbitrary element of the group U(2). It will be convenient to decompose 
into a product of a phase factor e~ ta £ U(l) and a special unitary matrix 
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ft' G SU(2): 

ft = e - ia n'. (3) 

This decomposition is ambiguous, as both the factors can be multiplied by — 1. 
However, the final results obtained with this decomposition will be free from this 
ambiguity, which makes the specific choice of the decomposition irrelevant. It is 
also worthwhile to note that in contrast to unitary transformations on the Hilbert 
space of quantum states, for which the overall phase factor is not physical, in the 
present case the overall phase factor e~ ia is physically meaningful, as it describes 
the phase of the fields which can in principle be measured with an external phase 
reference. 

The two-mode Hilbert space 7i describing our system has a convenient or- 
thonormal basis in the form of Fock states defined in general by 

\m H n v ) = i== — I vac), m, n = 0, 1, 2, . . . , (4) 

V mini 

where |vac) is the vacuum state of the system. The polarization transformation 
defined in Eq. ((J) does not change the total number of field excitations. It is 
therefore convenient to decompose the Hilbert space TC into a direct sum 

oo 

H = 0W« (5) 

1=0 

of finite-dimensional subspaces TiS 1 ' that contain exactly I excitations in both the 
modes. The subspace has dimension I + 1 and it is spanned by Fock states of 
the form: 

H {1) = Span{|m H (Z - m) v ) | m = 0, 1, . . . , I}. (6) 

Under the polarization transformation given in Eq. ijTJ). the state \mjj{l — m)v) is 
transformed according to: 

(4,r(4o'- m , 

\m H [l - m)v) = =Hvac) 
w ml [l — my. 

{Qhho'jj + n HV a\,) m (n VH a ! H + n vv a ] v ) l ~ m 

i — ^ . z 1 vac ) (t ) 

ym!(Z — m)\ 

If we now insert the decomposition of ft given in Eq. 0, this will produce an 
overall factor e~ tla times the same expression as in the second line of Eq. Q , but 
with the elements of fi replaced by those of ft'. It is easy to recognize in this 
expression the standard construction of irreducible representations of the group 
SU(2) using monomials |18| This yields the formula: 

l 

\m H (l - m) v ) - e~ Ua £ D l ^ n {n')\n H {l - n) v ) (8) 

n=0 
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where Dmn(fi') are the elements of (Z + 1) x (I + 1) matrices D l / 2 (fl') which form 
the irreducible (I + l)-dimensional representation of the group SU(2). Thus the 
unitary transformation U(fl) of an arbitrary quantum state in the Hilbert space 
H induced by the map given in Eq. (^Q) has the form: 

oo 

U(ft) = 0e- ,la ^ 2 (fi'). (9) 
1=0 

Let us note that although the decomposition ft = e~ ia ft' is defined up to — 1 
multiplying both the factors, the product e~ tla b l / 2 (ft') does not depend on the 
specific choice of the decomposition. This follows from the fact that the elements 
of the matrix D l / 2 {ft') are given by monomials of degree I constructed from the 
elements of ft'. Therefore each one of the substitutions e~ ta — > — e~ ta and ft' — ► 
— ft' will produce a factor (—1)', one multiplying e~ lla and another one multiplying 
D l / 2 (ft'), which will cancel each other. Consequently, the right hand side of Eq. © 
is defined unambiguously function of ft £ U(2). 



3. Collective depolarization 

We will now consider the scenario in which the polarization transformation ft 
is constant across ./V uses of the channel. The entire Hilbert space in this case 
is given by an TV-fold tensor product Tl® N of the two-mode space H analyzed in 
the previous section. The action of the collectively depolarizing channel on an 
arbitrary input quantum state \ip) 6 H.® N is given by: 

w ~ [u(fi)f N m- (io) 

Given the decomposition of U(ft) derived in Eq. ©, we can rewrite the tensor 
product [17(0)]®* as: 

[u(ft)f N = e - l(h+ - +lN ^b li / 2 {ft')® ...®b lN / 2 {ft') 

h,...,l N >0 
oo 

= ® e ~ iLa b h,2 {ft')® ...®b lN / 2 {ft') (ii) 

L=0 h,...,l N >0 
h+...+l N =L 

We see that the overall phase factor of the transformation ft enters the expression 
only with the total number of excitations L = l\ + . . . + Zjv contained in all the 
N uses of the channel. The A^-fold tensor product of the SU(2) representations 
b h / 2 {fl')®...®b l «/ 2 {fl') can in general be decomposed into direct sums by the 
iterative application of the formula [201 : 

b ji {fi) ® b j2 (ft') = D\&-to\(n') © b^- j2 \ +l {fi) © ... © b ji+j2 (fi'). (12) 
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This formula allows one to convert all the tensor products into direct sums of 
irreducible representations, with any of the representations allowed to appear a 
number of times. Therefore we anticipate that the inner sum in Eq. can be 
represented in the form: 

D h / 2 (n')®...®D l »/ 2 (n') = i Rj ®#(n') (13) 

h,...,i N >o j 

h+...+l N =L 

where Ik is the identity operator in a .fT-dimensional space C and the integer 
K 3 NL tells us how many times the representation j occurs in the decomposition 
of the iV-fold product [U(n)]® N in the sector of the entire Hilbert space Ti® N 
that contains exactly L excitations. The multipliticies K 1 NL will define the ca- 
pability of the system to protect quantum coherence against depolarization. By 
considering the parity of the indices l\, . . .In when applying Eq. (fT2|) to multiple 
tensor products D h l 2 {n') <g> . . . ® D iff / 2 (n'), it is easy to observe that the repre- 
sentations which will appear in the decomposition in Eq. (|13f) will be indexed with 
j = 0, 1, . . . , L/2 for even L and with j = 1/2, 3/2, . . . , L/2 for odd L. Assuming 
for the time being that the multiplicities K NL are known, we thus arrive at the 
following decomposition of the polarization transformation in the entire Hilbert 
space H m : 

oo L/2 

[U(m® N = ©e" !i ° i K , NL ®D j (tl'). (14) 

L=0 j=L/2-[L/2\ 

The above formula suggests the decomposition of the Hilbert space TC® N into 
sectors with a fixed number of excitations Hnl- 

CO 

^ = 0%! (15) 
L=0 

which in turn can be represented as isomorphic with the following structure: 

L/2 

n NL = c k nl®c 2j+1 . (16) 

j=L/2-[L/2] 

The components of the state vector belonging to separate subspaces Ti.NL are mul- 
tiplied by different phase factors e~ tLa . If the parameters of the transformation f2 
are unknown, this implies the loss of quantum coherence between different sectors 
Ti-NL- Within each sector, however, we have a number of subspaces for which the 
action of the collective depolarization operator is given by the identity 1 „j . This 
means that a quantum state which is encoded into a subspace isomorphic to one of 



[Author and title] 



6 



the subspaces C k nl in the structure defined in Eq. Q16JI remains intact even if the 
transformation fi is completely unknown. We note that the coherence between 
the subspaces C k nl with different values of j is destroyed as they are coupled to 
subsystems C 2 - 7 " 1 " 1 on which collective depolarization takes different forms. Con- 
sequently, the capacity of the system to protect quantum coherence is defined by 
the multiplicities K J NL , or to be precise, by the highest value over the permitted 
range of j. 



4. Recursion formula 



We will now demonstrate that the multiplicities K :I NL are related via a simple 
recursion formula. As we have seen in the preceding section, when the total number 
L of excitations is fixed, the overall phase of the polarization transformation fi is 
irrelevant. It is therefore sufficient to restrict ourselves to CI G SU(2) and consider 
the decomposition of [{/(fi)]®^ in the sector TLnl in the form: 



[I7(n)]« 



L/2 

© 

j=L/2-[L/2\ 



i s ®D>(n) 

J-^ AT T 



(17) 



The Hilbert space corresponding to N uses of the channel has an obvious decom- 
position as H m = H® {N - l) ®H. If we now consider the sector of Ti^ N containing 
exactly L excitations, it can be constructed from the subspaces of T^®^ -1 ) and TL 
according to: 



n NL = Hn-xjj ® n {L - L,) 



(18) 



L'=0 



where liS L is a subspace of the two-mode Hilbert space introduced in Eq. © 
that contains exactly L — L' excitations of the field. This construction implies that 
the operator \U(Cl)]® N _ can be represented as: 

[U{ft)f N 



-- [u(n)f( N -v 

L'=0 

We can now insert the decomposition of [[/(fi)]®^ -1 ) 

analogous t< 
This yields: 



,8^(^0/2(0) 



(19) 



using an expansion 



analogous to Eq. (fT7|) . assuming that the respective coefficients K J N _ 1 L , are known. 



Mm 



L>/2 



Hnl 



© © V ®^"(n) ®^ (L - L ' )/2 (n) (20) 

1 N-1,L' 1 



L'=0 \j'=L'/2-[L'/2\ 
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The tensor product '(«) ® l)( z '- i ')/ 2 (ft) appearing in the above expression can 
be decomposed using Eq. (|12|) into a direct sum according to: 

i'+(L-L')/2 

/3- ? ''(n)0/3( L ~ i ')/ 2 (n) = b\n) (21) 

i=|i'-(L-L')/2| 

Inserting Eq. (j^J) into Eq. (|2T7|) gives: 

L Z//2 j'+(L-L')/2 

[#(n)H =0 i^, ®#(n) (22) 

Wi L'=0j'=L'/2-LL'/2j j=|j'-(L-L')/2| iV ~ 1 ' L ' 

Comparing Eq. (|17|) with the above expression allows us to relate the multiplicities 
K J NL , describing decoherence-free subsystems for N uses of the channel, to the 
coefficients K~L_i y corresponding to the case when the number of channel uses is 

reduced by one. In order to find the explicit relation between K 3 NL and K J N _ 1 L i 
we need to change the order of summations over L' , j', and j in Eq. (|22jl such 
that the summation over j is the outermost and its limits do not depend on other 
summation variables. Then the inner summations over V and j' should yield a 
formula for K ] NL . 

The triple summation given in Eq. (|22|) can be visualized by plotting in a three- 
dimensional space the set of points (L',j',j) that are defined by the summation 
limits in Eq. (|22j) . The resulting grid is depicted in Fig. ^ First, we note that 
the conditions L' = 0, 1, . . . , L and j' = L' /2 — [L' /2\ , . . . , L' /2 define a triangular 
grid for the pairs {V ,f). This two-dimensional grid is shown with the help of dark 
grey points in Fig.^ For every pair (L',j') we have an allowed range of js defined 
by the limits of the third sum in Eq. ()22|) . The upper limit for j is always given 
by the plane specified by the equation 

j = f + {L— L')/2. (23) 

The lower limit for j is specified by one of the two conditions: 

j = (L-L')/2-f (24) 



or 

j = j'-(L-L')/2 (25) 

whichever gives a higher value of j. It is easy to verify that Eq. (|24|) is relevant when 
j' < (L—L')/2, whereas Eq. ()25|) gives the lower summation limit if j' > (L—L')/2. 
The three planes defined by Eqs. (|2.3|) -(|25 [) together with the fourth vertical plane 
specified by the condition j' = L' /2 form a tetrahedron which encloses the entire 
three-dimensional summation grid. 
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Fig. 1: The three-dimensional grid for summation variables (L',j',j) defined by 
the triple sum in Eq. I|22[). The allowed values of L' and j' form a triangular 
grid shown with dark grey points in the (L',j') plane marked as a light gray 
square. For every pair of V and j 1 the permitted values of j are shown as black 
points connected with a dotted line. The entire three-dimensional grid is contained 
within a tetrahedron with vertices located at (L,0, 0), (L/2,L/4, 0), (0,0, L/2), 
and (L, L/2, L/2). These points are obtained by calculating intersections of any 
three out of four planes given by Eqs. (j2^j) - (j25|) and j' = L' /2. 
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The structure depicted in Fig. ^ provides us with guidance on how to invert 
the order of summations in Eq. (|22|). in order to make the summation limits of 
j independent of other variables. Obviously, the summation over j will run from 
L/2 — \L/2\ to L/2 in integral steps. In order to find the limits for V and j 1 
for a fixed j, we need to consider an intersection of the tetrahedron depicted in 
Fig. ^ with a horizontal plane corresponding to that value of j. This procedure is 
illustrated in Fig. [2J The intersection of the plane of constant j with the planes 
specified in Eqs. (|23|) -(|25 |) gives three linear constraints on the values of V and 
j' shown in Fig. |2] with dashed lines. These constraints, together with j' = L'/2, 
define a rectangular region which, as it is easy to see, lies entirely within the trian- 
gular grid of pairs (L',f). Because of the geometry of this region, it is convenient 
to define two new variables: 

H = L'/2 + f, (26) 
v = L'/2-f. (27) 

The limits for /i are given by Eqs. (|24j) and (|25|) and have the explicit form L/2— j < 
[i < L/2 + j, whereas the values of v are bounded by j' = V /2 and Eq. (|23j) which 
combined together give < u < L/2 — j. Both the variables \x and v increase 
in unit steps, which follows from the form of the grid for the (L',j') variables. 
Consequently, the summation over V and j' for a specified j can be given by a 
double sum over fi = L/2 — j, L/2 — j + 1, . . . , L/2 + j and u = 0, 1, . . . , L/2 — j 
with the old variables expressed as L' = fj, + v and j' = (// — v)/2. 

Using the new summation variables fi and v gives us the following decomposi- 



tion for [U(fl)f 



with the inverted integration order: 



L/2 L/2+j L/2-j 

[u(nT N = i^, ®^(n) (28) 

j=L/2-[L/2\ n=L/2-j v=0 

The two inner sums yield of course the identity operator acting in a larger space, 
defined by: 

L/2+j L/2-j 

\J7 KA7 l^(M-")/2 = K (fl-v)/2 ■ (29) 

fi=L/2-j is=0 

Inserting this equation into Eq. ()28|) and comparing the resulting expression with 
Eq. 1)17(1 gives a recursion formula for K° NL in the form: 

L/2+j L/2-j 
fi=L/2-j u=Q 




Fig. 2: The points of the three-dimensional grid (L',j',j) that lie in the horizontal 
plane for a specified j have the first two coordinates (L',f) specified by the set of 
inequalities L'/2+f > L/2-j, L'/2+f < L/2+j, and L'/2—f < L/2-j. These 
three inequalities, together with the condition L'/2 — j' > 0, define a rectangular 
region, which can be conveniently parameterized with the help of two new integer 
variables n and v. 

which is the central result of this paper. In order to complete the recipe for calcu- 
lating the multiplicities Kir L , we need to specify their values for N = 1. This task 
however is straightforward. Let us recall that for N = 1 we have H N=1L = 
where the right-hand side has been defined in Eq. and that for a polarization 

transformation ft G SU(2) Eq. © implies that U(fl) = D L / 2 (£l). This 

means that for N = 1: 

K' N=ltL = 5 L , 2j (31) 

which provides the initial condition to calculate K J NL for an arbitrary N. 

The recursion formula derived in Eq. (j30|) can be illustrated with a diagram 
shown in Fig. |31 which provides a mnemonic recipe for carrying out calculations. 
In order to find K J NL we need to add K 3 N _ 1 L , for all pairs of (L',j') that lie within 



r 
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• \« 



(LJ) 



1 2 



L' 



Fig. 3: Graphic representation of the recursion formula derived in Eq. (|3U|) , The 
multiplicity K 1 NL is given by a sum of the multiplicities K J N _ l L , for all the pairs 
of the parameters (L',j f ) lying inside a rectangle marked as a light grey region, 
with vertices at (L,j), (L-2j,0), (L/2 - j, L/4 - j/2), and (L/2 + j, L/A + j /2). 

a rectangle which is rotated by 45° with respect to the axes of the coordinate 
system, and has two vertices located at (L,j) and (L — 2j, 0), with the remaining 
two located on the line j' = L' /2. Fig. |21 shows how to obtain this rectangle by 
starting from the point (L, j) and going in two orthogonal directions that form 45° 
with the axes, until reaching respectively the limits of j' = L' /2 and f = 0. 

5. Conclusions 

We have derived a recursion formula for the dimensions K J NL of decoherence- 
free subsystems in a bosonic channel experiencing collective depolarization de- 
scribed by the U(2) group. This model is relevant in quantum communication 
over single-mode optical fibers, for which collective depolarization is one of the 
dominant decoherence mechanisms. Although we have not been able to solve the 
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recursion formula and obtain a closed analytical formula for Kir L , it can be easily 
implemented in numerical calculations. 

The depolarization model considered here assumed that the phase relations 
between the consecutive uses of the channel are fixed. This requirement cannot 
be fulfilled if the communicating parties do not share a common phase reference 
to prepare and detect states. Then the phase factor e~ ia varies between the 
uses of the channel, and only the SU(2) transformation ft' remains constant. In 
order to implement decoherence-free encoding in such a scenario, the same fixed 
number of I excitations must be transmitted in a sequence of channel uses. The 
structure of decoherence-free subsystems for N uses of the channel is then given 
by a direct-sum decomposition of [D l / 2 (Tt')]® N , discussed in Alternatively, 
it is possible to devise schemes in which phases are self-referenced by employing 
multiport interferometers for state preparation and detection. Examples of such 
schemes have been described in |lf)(H2j. Mathematically, this approach consists in 
introducing phase dependencies between multiple uses of the channel, which then 
give rise to non-trivial decoherence-free subsystems. 
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